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Dissipation-free photon-photon interaction at the single 
photon level is studied in the context of cavity electromag- 
netically induced transparency (EIT). For a single multilevel 
atom exhibiting EIT in the strong cavity-coupling regime, the 
anharmonicity of the atom-cavity system has an upper bound 
determined by single atom-photon coupling strength. Photon 
blockade is inferred to occur for both single and multi-atom 
cases from the behaviour of transition rates between dressed 
states of the system. Numerical calculations of the second or- 
der coherence function indicate that photon antibunching in 
both single and two-atom cases are strong and comparable. 
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Enhancement of dissipation-free photon-photon inter- 
actions at the few photon level is a fundamental chal- 
lenge in quantum optics. Such interactions are necessary 
for single photon quantum control as in quantum logic 
gates. Typically, tuning close to atomic resonances si- 
multaneously increases both (Kerr) nonlinearity and lin- 
ear absorption, severely limiting the available nonlinear 
phase-shift. In Ref. Q, a new method that is based on 
electromagnetically induced transparency (EIT) [^ was 
proposed. In this scheme linear susceptibility vanishes 
on resonance, making it possible to obtain as much as 
8 orders-of-magnitude improvement in Kerr nonlinearity 
as compared to conventional schemes based on three-level 
atoms. Recently, experimental realization of this scheme 
with an enhancement approaching 10^ using a sodium 
Bose-Einstein-condensate (BEC) has been reported M. 
In light of these developments, it is natural to explore 
the limits on achievable photon-photon interaction while 
maintaining EIT. 

A single photon in a cavity can block the injection of 
a second photon due to a photon blockade effect [§,^. 
The original proposal for using a cavity-EIT medium to 
achieve photon blockade was based on an adiabatic elim- 
ination of the atomic degrees of freedom Q . Using a lin- 
earized analysis, Grangier et. al. M have shown that the 
adiabatic elimination carried out in Ref. ^ is not justified 
due to large dispersion of the EIT medium. The correct 
description of photon-photon interactions in the single 
photon limit requires a nonlinear analysis that treats the 
atom-photon interactions exactly: this is the principal 
goal of this Rapid Communication where the analysis 
is based on the exact energy eigenstates of the coupled 
atom-cavity system. Under conditions of strong coupling, 
the anharmonicity in the spectrum of the atom-cavity 
molecule can lead to nonclassical photon statistics such 
as photon antibunching or a photon blockade effect. The 



latter is determined by the probability of two photons 
absorbed sequentially into the system while the former is 
related to the conditional probability of photon emission 
given a single photon was already emitted. Both effects 
strongly depend on the anharmonicity of the atom-cavity 
energy eigenstates in the doubly excited (n = 2) mani- 
fold, or equivalently, how well can the coupled system, 
under EIT conditions, be described as an effective two- 
level system 0. Even though our analysis focuses on 
determining the level shifts and linewidths in the n = 2 
manifold, we complement this analysis with a calculation 
of the second-order coherence function of the transmitted 
field. 

We investigate the dependence of the photon-photon 
interaction on the atom number, dispersion, and the 
atom-cavity coupling strength. We concentrate on the 
situation where the cavity contains at most two photons, 
the atom-cavity system exhibits EIT at the cavity mode 
frequency, and a weak external driving field is resonant 
with the EIT transition so that the only state excited 
by the first photon (i.e. in the n = I manifold) is the 
cavity-EIT state. The crucial advantages of this EIT 
scheme are absence of one photon loss due to sponta- 
neous emission, arbitrarily reduced cavity decay rate for 
the n = 1 manifold cavity-EIT state, and a multiatom 
induced photon-photon interaction which, under condi- 
tions of non-resonant EIT, can be comparable to a single 
atom. In the context of quantum logic gates, atom-atom 
interactions between two multilevel atoms induced via 
a single cavity mode have been discussed previously |q| 
where the role of photons and atoms is reversed compared 
to our scheme. 

The definition of photon blockade suggests that its 
strength can be quantified by the transition rate VF1-+2 
from n = 1 manifold to n = 2 manifold, normalized 
to the n = — > n = 1 transition rate. Assuming a 
single relatively stable state in the n — 1 manifold, we 
can measure the strength of the photon blockade effect 
by -P = J2i^i^2i/Wo^i, which in turn suggests that 
one way to analyze the problem is to calculate transition 
rates between exact dressed-states. An experimentally 
more relevant quantity is the degree of photon antibunch- 
ing determined by g^{0) = {a^ d?)/{d)a)'^. Even though 
these quantities are physically different as described ear- 
lier, when the coupled system behaves like a single two- 
level system g^(0) <C 1 imphes photon blockade (P <C 1). 

Consider the n = 1 manifold of the coupled atom- 
cavity system where the cavity mode contains at most 
one photon. Since we assume throughout this paper that 
the EIT condition is satisfied, the frequencies of the clas- 



sical coupling field (wc) and the cavity mode {uJcav) are 
chosen to satisfy two-photon resonance tUcav — w^ — UJ21 = 
[^31 — 6 — ojc ^ UJ21 — (see Fig. 1). The EIT condition 
implies that there is no population in the upper state 
1 3) which is the only atomic state with a significant de- 
cay rate relevant to the n = 1 manifold. The effective 
interaction-picture Hamiltonian which couples N 4-level 
atoms with the single cavity mode and the coupling field 
with Rabi frequency 2Q is 



N 



H /Ti — —iKa'a + 2^ 
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+ 513(00-31 -I- a^d\^) + 524(00-42 -t- a^al^) 



(1) 



where 731 , 732 , 74 are atomic level spontaneous emission 
rates, r3 = (731 -I- 732)72 -I- J(5, r4 = 74/2 -t- iA, k is 
the cavity decay rate and A = W41 — UJ21 — cocav The 
externally measured photon statistics of the driven cavity 
are determined by the master equation p = —i[H/h, p] + 
£p \p) — a, p] -t- jp where j is the jump superoperator for 
Lindblad form and £p is the strength of the weak external 
drive field resonant with the cavity mode. The last term 
of Eqn.l can be replaced by Hperturb = —^X^^o,^i^22 
in the limit |A| ^ 524,1^4 where x = 524/^- 

In the case of a single atom in the cavity {N = 1), 
the n — 1 and n = 2 manifolds (in the perturbative 
limit) have 3 eigenstates. In the limit r3 = 0, eigenvalues 
Tiiujcav + Cj) for the n — 1 manifold are approximately 
given by 
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-i-±Qy^l + a-{K/2n)^ (2) 



where a = gf^N/il'^ and i?e(eg^^) is independent of 
{gi3,fl) and is on cavity resonance. The eigenstate cor- 
responding to Eq"^ 



\^r') 



VT 



gi3 



-0-21 



|o) 



(3) 



contains no component from the upper atomic state |3) 
and in the limit a 3> 1 is predominantly atomic in charac- 
ter: this is the cavity-EIT trapping state. The transition 
rate induced by a weak probe from the ground state into 
Cq^^ is independent of a. In the n = 2 manifold, in the 
limit 924 = 0, one of the eigenstates is 
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whose linewidth due to cavity decay has a weak a de- 
pendence of k(4 -|- a)/{2 + a). For 524 7^ 0, this energy 
eigenstate experiences an ac-Stark shift. Once the EIT 
supporting state \(f>o^^) is excited, the state closest to 
two-photon resonance for (5 = (|0o='^)) ^^^ ^ perturba- 
tive level shift due to 924 of 
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allowing photon-blockade and photon antibunching to oc- 
cur, provided x > i^- Level shift is largest in the limit 
a ^ 1, which can be achieved by reducing f2. One should 
note that additional constraints on f2 exist to maintain 
EIT but are weak for cold alkali atoms |3|. 

The results outlined above only consider the regime 
where level |4) can be adiabatically eliminated. In this 
perturbative limit, we found that the maximum value 
of the photon-photon interaction (as measured by the 
anharmonicity or level shift in analogy with an ideal 
nonlinear cavity) is independent of a and is given by 
\Xraax,pert\ ~ |x|- This result suggests that the upper 
bound for the non-perturbative photon-photon interac- 
tion will be given by 324 in the single-atom case. Consider 
energies h{2ujcav + e) of the n = 2 manifold eigenstates. 
In a nonperturbative single atom regime, there is no ev- 
idence that detunings increase the anharmonicity. The 
limiting case of vanishing detunings and losses gives: 



^ = ±^^±)J^-29h9h 



(6) 



where G^ = g|4 -I- fi^ -I- 2gi^. For the resonant case 
(A = (5 = 0) with g — 1713 = (724, the eigenvectors can 
be found exactly and eigenstates whose energy level is 
closest to 2TiLJcav have level shifts which asymptotically 
approach g from below as g increases. We believe that 
this result strongly suggests that the maximum possible 
anharmonicity achievable in cavity-EIT system is always 
bounded above by the single atom-cavity coupling co- 
efficient (724 • When 324 ~ V2513 the non-perturbative 
nature shows itself through level anti-crossings. There is 
then a second upper bound given by v25i3- 

In order to relate these ideas to well known photon 
correlation measurements we calculate 5^(t) for experi- 
mentally achievable single atom parameters (Fig. 2) [^. 
For weak cw driving such that the average cavity photon 
number is much less than unity, we expect g^(0) ^ 1 
for sufficiently large energy shifts in the n — 2 mani- 
fold. If the coupled system under weak excitation be- 
haves ideally like a single two-level system then no am- 
plitude for photon modes with n > 1 will contribute to 
.9^(0). One can understand this simply as restating the 
fact that photon blockade implies that ideally, photons 
transmitted out of the cavity should only correspond to 
cavity decay of |0q^^). Components from other transi- 
tions can be strongly suppressed by large splittings of 
the dressed states. The main contribution for nonzero 
3^(0) would then come from |0q=^). However, the val- 
ues we find for g^{Q) are as small as 2x10^^ for pa- 
rameters from recent experiments ||lO(| on cavity-QED 
{g = 7.5k, 731 = 0.325k). Such experiments might be 
enhanced using EIT. For Vt = 2.5731 so that a ~ 10^ one 
would expect a reduction of the cavity-EIT decay rate 
from their measured cavity decay rate of 16 MHz to 0.19 



MHz. Furthermore, future improvements in cavity and 
atom trap design will lead to cavity-QED experiments 
limited by the atomic spontaneous emission rate; cavity- 
EIT schemes will be especially favourable in such a case. 
Clearly the ultra- low values of 17^(0) can not be explained 
by the level shifts alone. Since the weak probe is tuned 
exactly between the two dressed states in the n — 2 man- 
ifold quantum interference might be responsible. 

We next discuss how in a high dispersion regime the 
usual linear scaling of the nonlinearity with the number 
of atoms does not apply. The n = 1 manifold, where a 
single photon excitation is shared among the N atoms, 
contains three energy levels h{ujcav + ej) whose energy 
shifts and widths, neglecting spontaneous emission, are 
approximately given by Eqn.(2) provided we assume that 
all atoms couple identically to the cavity mode. Note that 
the eigenvalue Cq^^ has a width which decreases with 
the number of atoms N |1^ while the splitting of the 
other two dressed states increases with N. The eigenstate 
corresponding to e^^^ is 
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When A^ > 1, the n = 2 manifold contains six energy 
eigenstates (in the perturbative limit). For 524 = 0, one 
of these eigenstates 
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has energy h2ujcav, implying that the energy level dia- 
gram, as probed by the weak external field, is harmonic. 
The state jf/)"^^) corresponds to two cavity-EIT trapping 
state excitations and has a linewidth which decreases 
with 1 -|- a. For (724 7^ 0, this energy eigenstate expe- 
riences an ac-Stark shift. When there are no degenerate 
eigenstates the level shift of \(f>"^'^) caused by the pertur- 
bation Hperturb is 
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The atom-cavity molecule therefore has an anharmonic 
response to the drive field, provided that |e""^ | > 2k/(1 + 
a). We note that the (nonlinear) splitting of the eigen- 
state corresponding to e"=2 initially increases with the 
number of atoms N in the cavity as would be expected 
from a traditional nonlinear optical system. However, as 
dispersion becomes important (a ~ 1), this increase sat- 
urates. In the highly dispersive limit given by a ^ 1, 
the splitting of the eigenenergy e""^ and its width de- 
creases with increasing N, under the assumption that all 
TV atoms have the maximum possible interaction. 

As we have already argued, the photon-photon interac- 
tion strength is determined by the energy level splittings 
and linewidths in the n = 2 manifold, as well as the 



transition matrix elements. In the special case of doubly 
resonant FIT {S = 0) with iV > 1, one of the eigenstates 
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has energy h2ujcav and a linewidth that is independent 
of a. This state remains unshifted by 324 and has a 
nonzero transition rate from It/)"^^) that decreases with 
a. Degenerate perturbation theory shows that the ac- 
tual shifted eigenstate is then a linear combination of 
|V'"=^) and |(/)"=2) with level shift 3e"=V2. The tran- 
sition rate from the n = 1 manifold FIT state to this 
shifted eigenstate increases with the number of atoms iV, 
for a > 1, X > Oii^- Increasing the number of atoms there- 
fore leads to a reduction in the magnitude of the photon- 
photon interaction. The presence of these two states im- 
plies that photon blockade is small for this special case of 
i5 = 0, where one would naively expect the nonlinearities 
to be the strongest. It is also interesting to note that the 
detrimental effect of the degenerate state |V'"^^) that in- 
hibits nonlinearity implies that the adiabatic elimination 
of atomic degrees of freedom is not justified even in the 
low dispersion limit where it is normally assumed to be 
valid. 

When |e"=2|,K < |(5| < ^/Wgl^'+V^ , the degeneracy 
in the n = 2 manifold is removed and the dominant con- 
tribution for transition from n = 1 to n = 2 manifolds is 
determined by |0"^^). In this case, P becomes indepen- 
dent of a for a ^ 1 and approaches the single-atom re- 
sult. Photon blockade experiments using the cavity-EIT 
system with A^ ::^ 1 can therefore achieve a reduction 
in two-photon absorption by the atom-|-cavity molecule 
comparable to that of a single atom-|-cavity. However, 
the multiatom case is considerably more sensitive to other 
effects not included in the dressed-state calculations: for 
example, effective Rabi frequency associated with the ex- 
ternal drive has to be much smaller than e"=^, which 
places a constaint on £p. Also, the spatial dependence 
of the atoms, their finite velocities or atom-atom inter- 
actions may become more important. A recent analysis 
of the spatial dependence of atomic beams ||l^ for the 
Jaynes-Cummings model(JCM) where energy levels are 

proportional to Gjcm = 



E^i(ff^/5' 



probability distribution P{G 



JCM) 



)2, shows the 
is well approximated 



by a Gaussian with mean ^ N^ff , given by the atomic 
density and cavity parameters, when N^ff ^ 1. In the 
limit Neff ^ 1 and a ^ 1, the variation in e"=2 jgyel 
splitting and linewidth over the width of the distribution 
P{Gjcm) both decrease like 1/a^. The cases N^ff ~ 1 
and a ~ 1 are more involved and will be reported else- 
where. 

The structure of the eigenstates and energy eigenval- 
ues is the same for any A^ > 1 and so we have calculated 
(7^(0) for N — 2 versus detuning 5 of the cavity mode 



from UJ31 while maintaining EIT using the couphng con- 
stants and decay rates as in the single atom case (see 
Fig. 3). The energy level splittings in the n — I man- 
ifold constrains the detuning so that 6 < ^Jgl^N -\- fi^; 
we find that large 5 also results in suppression of energy 
level splittings in the n = 2 manifold. As can be seen in 
Fig. 3 values for g^{Q) using two atoms varies by three 
orders-of-magnitude over a detuning 0{g); we find that 
for such detunings 5^(0) ^ 1 and its value is of the same 
order-of-magnitude as for a single atom for the same driv- 
ing strength. In contrast, when (5 = 0, there is practically 
no photon antibunching for the two-atom case - a result 
we anticipated based on the appearance of the unshifted 
eigenstate |V'""^)- Furthermore, for \6\ = g, £p = 0.1k, 
fi in the range 2k ~ 10k we find g^(Q) — 10~^. As long 
as the driving satisfies £p < k/{1 + a) then large an- 
tibunching should be robust. These values for g^{0) are 
about two-orders-of-magnitude smaller than for similarly 
prepared two-level atoms in a cavity. 

In summary, we have shown that photon-photon inter- 
actions using cavity-EIT is experimentally feasible, and 
has distinct advantages over two-level atoms. We used 
the eigenstates of the perturbed many 3-level atom sys- 
tem to give the A'' atom dependence of the energy level 
shift. Using a heuristic argument based on transition 
rates between dressed states we inferred that multiatom 
photon blockade occurs provided non-resonant EIT is 
used. We have shown explicitly that photon antibunch- 
ing for a multiatom case (two atoms) can be comparable 
to that a single atom. We also showed that in the single- 
atom case, anharmonicity is bounded above by the single 
atom-photon coupling strength 524 • This bound is not 
reached if a/2(7i3 < (724 • 
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FIG. 1. (a) Internal atomic degrees of freedom of a 
four-level atom coupled to a single cavity mode with frequency 
iJcav and an external laser field characterized by the Rabi fre- 
quency 2f2. (b) The energy levels of the coupled atom-cavity 
system. 



FIG. 2. Second-order photon correlation (/^(r) demonstrat- 
ing photon anti-bunching. The coupling constant and de- 
cay rates correspond to current single atom experiments [10] 
of 513 = 524 = 7.5k, 74 = 731 = 732 = 0.325k and 



6 = A = 0,n = 2.5'y3i,£p = 0.1k. 



FIG. 3. Two atom case (A'^ — 2) showing the variation of 
g^{0) with the detuning S = UJ31 — ujcav for the cavity-EIT 
giant Kerr nonlinearity. The coupling constant and decay 
rates are the same as in Fig. 2 while the detunings and drive 
strengths are A = 0, £^p = 0.01k. 
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